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Confinement and chiral symmetry breaking are the two main facts in low-energy 
QCD and show up explicitly in the spectrum of quark bound states. Confinement im- 
plies that the nature of the quark dynamics has to change strongly in dependence on 
the scale of energy or distance associated with the interaction. Let us consider the 
interaction between a quark and an antiquark. At short quark-antiquark distance r, 
the interaction is dominated by the one gluon exchange (dipole configuration) and is 
obtained by a perturbative expansion in a s . With increasing distance, the interaction 
becomes linear in r (flux tube configuration). The formation of a flux tube of constant 
energy density o ~ lGeV/fm (the so called string tension) confining the quark and the 
antiquark to remain bound, is the most striking feature of confinement and is inher- 
ently nonperturbative. We have hints about this flux tube from the phenomenological 



spectrum and the Regge trajectories. Phenomenological potential models include in- 
teractions rising with the quark- ant iquark distance. Dual models explain this feature 
as a chromoelectric flux tube formation due to a dual Meissner effect, the QCD vac- 
uum behaving like a t'Hooft-Mandelstam superconductor. Confinement is then due to 
monopoles condensation. 

Lattice studies of the electric and magnetic fields strength distributions around 
a static quark- ant iquark pair yield direct evidence and a detailed picture of the string 
formation (see Figfj]). Lattice data provide us with: a quark-antiquark static potential 
linear for r > 0.2 fm; an interquark flux tube structure with a transversal rms width 
ranging between 0.5 fm and 0.75 fm; positive tests of monopole condensations; Abelian 
projected electric fields and currents distributions in agreement with the dual Ginzburg 
Landau predictionsBS 

The lattice remains the most suitable tool to obtain nonperturbative information. 
However, most of these appear as independent and uncorrelated measurements. On the 
other side we have many phenomenological models implementing confinement. Finally, 
high order perturbative calculations are also available. In order to gain some insight in 
the mechanism of confinement and flux tube formation we need a mo del- independent, 
exact (in the sense that it contains perturbative and nonperturbative physics and is 
valid at any scale), gauge- invariant and systematic approach to the quark interaction. 
It helps to have a formulation which is closely related to the lattice, allows analytic 
calculations and is physically transparent. 

In Sees. 2 and 3 we show that this program is realizable in the heavy quark case 
since a small expansion parameter (the quark velocity v) still exists. The result obtained 
is model independent and suitable for lattice evaluation. An analytic evaluation is 
possible once a QCD vacuum model is considered. Therefore it is possible to test models 
of confinement in an unambiguous and direct relation with the phenomenological and 
lattice data. 

We know that in the low energy region new degrees of freedom become relevant. 
The flux tube formation hints to the fact that dynamical gluon effects are vital when 
attempting to incorporate nonperturbative interaction, indeed the gluonic degrees of 
freedom are condensed into stringlike flux tube. Hence, the flux tube should determine 
the nonperturbative bound state quark dynamics. On one side, this means that also 
excited states of the interquark string has to be considered. This leads to the consider- 
ation of hybrids states and was discussed by Michaelu at this Conference. On the other 
side, this means that even when considering the ground state of the interquark glue, 
the flux tube energy should be taken into account properly. Moreover, to understand 
the transition region, we should be able to control precisely how the collective behavior 
of gluons modifies and alters single gluon exchange. 

In Sees. 4, 5 and 6 we discuss the heavy quark dynamics in QCD vacuum models 
that reproduce the flux tube structure, in comparison with the lattice data and with 
the "l/Q 4 " models extensively used in the literature. We find that all the 'flux tube 
models' predict the same form for the heavy quark interaction also in the intermediate 
region. We show that the Gaussian dominance of the Wilson loop average has to be 
related with the Abelian dominance of infrared QCD establishing a connection with 
the dual superconductor mechanism. 



In Sees. 7 and 8 we extend our discussion to bound systems involving at least one 
light quark. Here the spontaneous chiral symmetry symmetry breaking emerges in an 
interplay with confinement. 

2. Gauge-Invariant approach to quark bound states 

We consider a gauge invariant quark-antiquark singlet state[] 



>>J = ^i(x)U ik (x,y,C)4(y)\0) (1) 



where i,j,k,l = 1, ...3 are color indices and the Schwinger string line U has been 
inserted to ensure gauge-invariance 

U(x,y;C) = Pexp{ig f A,(z)dz»}. (2) 

Jy,G 

At a time t — 0, a quark and an antiquark are created, interact while propagating for a 
time time t — T at which they are annihilated. Then, the quark-antiquark interaction is 
contained in the gauge- invariant four-point Green function (xj = (xj,T), yj = (y^, 0)) 

G(x 1 ,x 2 ,y 1 ,y 2 ) = i J VATrS(x 2 ,y 2 ; A)U(y 2 , yi )S(y u x i; A)U{x x , x 2 )e iS ^ (3) 

where we neglected the fermionic determinant and the annihilation termsfj], Tr is the 
color trace and S(y,x;A) is the quark propagator in the external field A. To find a 
closed expression for the Green function in Eq.(|3|) is quite a formidable task. However, 
the situation becomes simpler if we consider the case of heavy quarks. In the static 
limit m — ► oo the equation for the static quark propagator 

(i>y D -m)S (x,y;A) = S\x - y) (4) 
is solvable in closed formiS and we findi ii 

G ( T )"!^V (xi _ yi)5 s (x2 _ y 2 ) e -^^ T (W(T )) (5) 

where 

(W(T )) = VATrPexp{ig f dz»Ap(z)}ew{igSYu} (6) 

is the static Wilson loop average, the loop T being a rectangle. From the Feynman-Kac 
formula and Eq. (|5|) we get 

V (r)= lim Un(W(V )). (7) 

1 — >oo ]_ 



* Actually this state is not an eigenstate of the Hamiltonian. It serves as a trial state to extract the 
energy of the lowest eigenstate having a non- vanishing projection on \<f>). This is equivalent to say 
that the contribution of the string is negligible in the limit T — > oo, which will become more apparent 
when dealing with the Wilson loop. 

t Since we neglect the fermionic determinant we work in the quenched approximation. In the casejpf 
heavy quarks however we can introduce into the effective action the expansion of the determinant B 



The static quark-antiquark potential, if exists, is given in terms of the static Wilson loop 
average. Notice that this is an exact expression that contains both the perturbative and 
the nonperturbative dynamics. Indeed, no expansion in the coupling constant has been 
performed. Here and in the following sections we implicitly assume the existence of the 
potential and we postpone this issue to Sec 6. The spin- dependent 1/m 2 corrections 
to the potential were calculated inEM considering the first correction to Eq. @ and 
evaluating it on the zero order exact solution. The result is given in terms of vacuum 
expectation value of (color) electric and magnetic field insertions in the static Wilson 
loop and controls the fine and hyperfine separation in quarkonia. 
However there were some residual difficulties: 

• The kinetic energy terms was not considered^ lacking in consistency with the 
virial theorem. 

• The spin-dependent potentials appeared as an expansion in 1/m 2 completely 
missing the terms in lnm obtained at one loop in the S matrix perturbative 
calculation of the interaction.^ 

• The calculation of the nonperturbative contributions was based on several as- 
sumptions and there was no clear and unambiguous procedure to calculate the 
nonperturbative behavior of the v.e.v. of the field strength insertions in the static 
Wilson loop. 

• It was not clear when a potential description was holding. 

The solution is to disentangle the different scales of the bound system (m, p = mv, 
E = mv 2 , v being the heavy quark velocity) and to establish a systematic , unambiguous 
and exact expansion procedure. This will be reported in the next section. 

3. Analytic closed expression for the heavy quark interaction 

In the heavy quark systems, the existence of an expansion parameter (the inverse 
of the mass m in the Lagrangian and the velocity v of the quark as a dynamical defined 
power counting parameter) makes possible to establish a systematic expansion proce- 
dure. The tool is provided by NRQCD.E3 This is an effective theory equivalent to QCD 
and obtained from QCD by integrating out the hard energy scale m. The Lagrangian 
comes from the original QCD Lagrangian via a Foldy-Wouthuysen transformation. 
The ultraviolet regime of QCD (at energy scale m) is perturbatively encoded order by 
order in the coupling constant a s in the matching coefficients which appear in front 
of the new operators of the effective theory. This ensures the equivalence between the 
effective theory and the original one at a given order in 1/m and a s . At order 1/m 2 
the NRQCD Lagrangiandescribing a bound state between a quark of mass m% and an 
antiquark of mass m 2 iai§EJ 

D 2 m D 4 m a B m D E-E D 
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* The static quark propagator is the zero order solution. This is similar to the HQET approach. 
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+ ic s g 5 ipi + antiquark terms (1 2) H WiW^Wk^i 

8mf J m\m2 

+ ^^1t> 2 Q^i + -^4T a <Tii>2il>lT a afa. (8) 
mim 2 mim 2 m!m 2 

This is the relevant Lagrangian in order to calculate the bound state energies up to 
order 0(v 4 ). The coefficients evaluated at a matching scale 11 for a 

particle of mass rrij. The matching coefficients are 1 or at the tree level. At one loop 
they contain terms in In rrij which establishEJ the agreement of the present calculation 
with the result for the perturbative interaction obtained at one loop in the S matrix 
formalism.0 For simplicity in the following we consider the tree level values |§| 

With Eq. (H) we have reduced the heavy quark dynamics to a nonrelativistic 
Schrodinger theory for the heavy quark (antiquark) with Pauli propagator K coupled 
to the usual relativistic field theory for light quark (omitted here, quenched appr.) and 
gluons. This result is especially useful for calculation on coarse la±ti£e.llHl However, here 
we are looking for analytic results. The strategy is the following .ll§0 At the 0(v 2 ) the 
Pauli propagator satisfies the equation 
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m J + KzriPj ~ 9&) 2 - ^3 (Pj - gA) 4 + gA° (9) 
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and therefore admits a path integral representation 

K 3 (x,y\A) = J ^/,.i» ( ;T//'""'^ "' (10) 

where H is the Hamiltonian appearing in @. Starting, as in Sec. 2, from the qq Green 
function (||]), we can substitute the quark Dirac propagator with the Pauli propagator 
and use Eq. ([H]) to obtain the closed path integral representation 



G {T) — K{T) = j P[z 1>Pl ] J D[z 2 ,p 2 ]exp 
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x ^TrT s Pexp{^jf dz? A^z) + (11) 

% I i - ^ lkr ^ -^ dvf M) 

with (/[-A]) = \ J VA f[A]exp{iSyM} , T s spin ordering and P path ordering. Notice 
that now the kinetic energy of the quark is explicitly considered and the contour integral 
in A^ and F^ u is extended to the distorted Wilson loop in Fig. ||. Indeed, due to the 
presence of the path integral sum, we are considering any possible trajectory for the 
quark (antiquark) at variance with the static path of the previous section. 

§For a discussion of_the matching coefficients, the 4-fermions interaction and the power in v of the 
operators in M) seeB 




Fig. 2. Distorted Wilson loop. 



Then, the potential exists if we can make the following identification in 
I. 
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From Eqs.flTTl) and (|l~2| ) we findS 
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where {(f(A))) = (f(A)W(T))/ (W(T)) is the vacuum expectation value in presence of 
a quark-antiquark pair (i.e. the Wilson loop). Eq. ([13]) is the quark-antiquark potential 
at order v 2 of the systematic expansion in v. No expansion in the coupling constant has 
been performed. The result is physically transparent. In Eq. ( |i~3"D the first term contains 
the static and the velocity-dependent potential, the second is the magnetic interaction, 
the third is the Thomas precession, the fourth is the Darwin term and the last one 
is the spin-spin interaction. All the dynamics is contained in the Wilson loop and in 
the v.e.v. of field strength in presence of the Wilson loop, hence is pure gluodynamics. 
Had we worked in QED, we would obtain the same result as in ([13]). Notice that the 
coupling underlying (|T3|) is the vectorial coupling of the Dirac equation. The difference 
between QED and QCD is contained in the behavior of the v.e.v. of field strengths 
and Wilson loop. The non linear dynamics of QCD determines the nonperturbative 
behavior of these v.e.v and accounts for the difference with QED. 

Varying the quark path Zj(t) — > Zj(t) + 5zj(t), the Mandelstam formula can be 
obtained 



(f&FrWF^te))) - ((F, v { Zl )))((F Xp {z 2 ) 



5S x p{z 2 



with SS^ u (zj) = dzj5zj — dzj5zj. We conclude that to obtain the complete quark- 
antiquark order 0(v 2 ) interaction (quenched) no other assumptions are needed than 
the behavior of (W(r)): given (W(T)) everything is analytically calculable. On the 
other side, expanding the average on the distorted Wilson loop T in terms of the static 
Wilson loop T we get expressions for the potentials suitable for lattice evaluation.lia 
Then, the way to validation of analytic models of the QCD vacuum via lattice data and 
phenomenological data is open. The aim is to obtain information as much as possible 
model-independent on the nonperturbative dynamics. 

3.1. General form of the potential 

At the order 0(v 2 ) and with tree level matching coefficients the quark-antiquark 
interaction reads 

V = V + V SD + VVd (14) 
with the spin-dependent interaction 

= sfe? + 4) + A mr) + Ur)] pssr " ^sr) ; W(r) + 2v:(r)] (15) 



1 (t c t C \^'( r ) , 1 /Si-rS 2 -r Si • S 2 \ . , , Si • S 2 . , , 
+ Li • S 2 - L 2 • Si + 5 q — v w) + q V ± r 

m\m,2 r mim 2 \ r z 6 J 6mim 2 

where Ey = r x p^, and the velocity-dependent interaction is given by (| | = Weyl 
ordering) , 



Wd = {pi ■ PaVb(r)} + 

mim 2 




(16) 



The dynamics is contained in the Vi(r) functions. These can be evaluated analytically 
via functional derivatives or expansion of the_W.ilson loop (once its nonperturbative 
behavior is assigned in a vacuum model), seeu'u or can be evaluated numerically—on 
the lattice from v.e.v. of field strength insertions in the static Wilson loop, see.liS'EJ 
The spin-dependent potentials agree wj±hthe Eichten-Feinberg result .0 The velocity 
dependent potentials were obtained inEJ'ta and have to be taken into account to be 
consistent at the order v 2 of the systematic expansion. The lattice results are presented 
in Sec. 5 while in the next Section we summarize the analytic results obtained in various 
models of the QCD vacuum. The common thread of these models is reproducing the 
interquark flux tube structure. 

4. Infrared dynamics and flux tube structure 



The expression ( |13[ ) for the heavy quark interaction is exact at the order of the 
expansion in v considered. It is gauge-invariant and this allows us to insert approxi- 
mations, i.e. vacuum models to get the nonperturbative behavior of the Wilson loop 
average. We retain the relevant configuration in the nonperturbative regime, which 
turns out to be a flux-tube like configuration. 



4-1. A simple model calculation 



The simplest model consist in assuming that the short and the long range con- 
tributions to the Wilson loop factorize and evaluating the last one using an area law 
(as suggested by the strong coupling expansion) 

\og(W(T)) = \og[(W(r)) SR -(W(r)) LR ] = -^ £dx* fdtfD^x-y^aS^ + ^V. 

(17) 

D^v is the free gluon propagator, a is the string tension, S'min is the minimal area 
enclosed in the distorted Wilson loop contour (cf. Fig. |^), C is a constant and V is the 
perimeter of the loop. The generaLexpression for the minimal surface with contour V 
is the Nambu-Goto action. In ref.ji3 it is shown that, at the order v 2 of the interaction, 
the minimal area is equivalent to the area spanned by straight lines connecting pc 
at equal time on the quark and antiquark trajectories. For any details we refer to.E 
Here, we only present the final result of this minimal area law model (MAL) 

\r { \ 4 «s . < n \r ( \ 8as ° \r( \ 2a * a t/ t \ a C 

o(r) = "77 + or + b(r) = "97 ~ 9 r c(r) = ~~37 ~ 6 r d(r) = "9 r ~ 7 

V e (r) = ~r V[{r) = -o V>(r) = ^ V 3 (r) = 4^ V 4 {r) = f 7ra s 5 3 (r) (18) 

and AV^ = 0. The one gluon exchange term (proportional to a B ) coincides with the 
usual Breit-Fermi potential. Let us focus on the nonperturbative part (in a and C). The 
V1 — V4 spin potentials (sd) agree with the electric confinement calculation of Eichten 
and FeinbergQ (taking into account the Grpnies correction!!!]) . The V a — V e velocity 
potentials (vd) were first calculated in MAL.EJlla To make clear the physical content of 
the result, let us write the nonperturbative contributions in the center of mass system 
for equal masses. The spin-dependent part is 

and corresponds to pure Thomas precession, i.e. the magnetic contribution is zero. The 
velocity dependent part is 

V™ = ~^K- ( 20 ) 
6 m z r 

and it is proportional to the flux tube angular momentum. Indeed, the same vd cor- 
rection is obtained in the relativistic flux tube modelEiJ where the energy of the in- 

terquark flux tube is explicitly added to the HamiltonianQ. Then, it is clear that our 
gauge-invariant Wilson loop approach has automatically included the energy of the 

flux tubejjj]. 

^In the relativistic flux tube model the Lagrangian is composed by the relativistic energy of the 
quarks plus the relativistic energy of the flux tube considered as-a mechanical tube with constant 
energy density and velocity transversal to the interquark line, see.Ell 

II However, we are not using any 'mechanical model.' 



The result is even more interesting since it cannot be obtained within the semirela- 
tivistic reduction of a pure convolution (i.e. depending only on the momentum transfer 
Q = pi — p[) Bethe-Salpeter kernel, with any mixture of vector and scalar components. 
In general, a scalar 1/Q A BS kernel is used in order to reproduce ( fL9f) , however the 
velocity dependent correction differ from (p0|), see Sec. 5. On the other hand, the re- 
sult (|19|)-(|20D is physically transparent. Imagine a quark-antiquark pair connected by 
a chromoelectric flux tube. The magnetic field is zero in the comoving frame and then 
the spin-interaction is purely Thomas precession; the electric tube is moving with a 
transversal velocity and then its energy originates the vd term. This is the Buchmiiller 
picture. We emphasize that the 'scalar-like' character of the spin splitting we obtained 
is dynamically generated through the collective nature of the gluonic degrees of free- 
dom. It is interesting that inc3 this same conclusion has been obtained working in the 
framework of the Coulomb gauge Hamiltonian and diagonalizing a sector of the Fock 
space. Again, the scalar interaction is effectively generated. However, the Wilson loop 
approach appears to be simpler and more powerful (at least for heavy quarks); the 
intermediates states that are explicitly considered in the Fock Hamiltonian approach 
here are summed up into the functional formalism. Expressing the potential in terms 
of the Wilson loop we keep the relevant degrees of freedom. Our special 'trial func- 
tion' constructed with the string interquark operator has selected out the flux tube 
configuration. 

The MAL model is in many respects too rough (e.g. the {{F^ u F pa )) v.e.v. turns 
out to be zero and the flux tube turns out to be infinitely thin) and it may play some 
role only in the limit of very large interquark distances where an effective relativistic 
string model description is expected to hold. In the next Sections we discuss more 
sophisticated assumptions. 

4-2. Gaussian dominance in gluodynamics versus Abelian dominance in infrared QCD 

We need models of the QCD vacuum which provide us with the nonperturbative 
behavior of the Wilson loop average. To this aim we want to exploit all the available 
lattice information on the mechanism of confinement and the measurements of Wilson 
loop andJIeld strength v.e. v.. About the mechanism of confinement, at this Conference, 
evidenceofl was presented that the QCD infrared dynamics is well approximated by a 
(dual) Abelian Higgs model. Indeed, in the Abelian projection the QCD gluodynamics 
is reduced to Abelian fields, Abelian monopoles and charged matter fields degrees of 
freedom. Lattice simulations show Abelian dominance and monopole dominance in the 
long range features of QCD and condensation of monopoles in the confined phase. For 
details we refer to.30 The measured electric fields and magnetic currents in the presence 
of static quark sources are consistent with dual Ginzburg-Landau type of equations. 
The penetration length A = 1/M (M = dual gluon mass) and the correlation length 
£ = 1/Mfj, {Mfj, = Higgs mass) are measured. It is found M ~ and therefore the 
QCD vacuum behaves as a dual superconductor on the border between type I and type 
II. Flux tube solutions (of the type of Abrikosov-Nielsen-Olesen vortices) exist and the 
structure of these solutions is controlled by the penetration and the correlation length. 
However these are lattice measurements and it is not clear e.g. the relation of A and £ 
with the parameters of QCD. 



On the other side, let us consider the v.e.v. of the Wilson loop. It pays to expand 
this average in terms of fields strength expectation values, by using the non-Abelian 
Stokes theorems 



(W(T)} = (exp{ig I dz^{z)}) = (Pexptg f ^(1)^(0,1)^(1)^(1,0)) = (21) 
ex P{E d <Wl) • • -^n(n)(^(0, 1)1^(1)17(1, 0) • • .F^(n)l7(n, 0)) cum } 



n=0 



where i ee Xj, S(T) denotes a surface with contour T and (. . .) cum stands for the 
cumulant average.Bj In principle, the v.e.v. of the Wilson loop is given by the sum of 
all the cumulants. However, in a recent lattice investigation,^ evidence of the Gaussian 
dominance in the cumulant expansion of quasi-static Wilson loop average was found. 
Therefore, it follows that 

(W(T)} ~ exp{-i J s r dS^O) dS pa (x) (g 2 U(0, x)F^(x)U(x, 0)F pCT (0))} (22) 

is a good approximation. This is the basic assumption in the (Gaussian) stochastic 
vacuum modeled and it was phenomenologically confirmed by calculation in high energy 
scatteringc3 and quarkonia. Then, the heavy quark interaction is determined by the 
two-point field strength correlator 

g 2 (U(0,x)F pu (x)U(x,0)F Xp (0)) = {{5, x 5 vp - 5 pp 5 uX )(D(x 2 ) + D 1 (x 2 )) 

+{x il x x 5 vp - x p x p 5 uX + x u x p 5 pX - x^x^^-^D^x 2 )^. (23) 



In (23) the Lorentz decomposition is general and the dynamics is contained in the form 
factors D and D±. The function D is responsible for area law and confinement (indeed 
in QED, Hue to the Bianchi identity, we have D — 0). For D and Di the lattice cal- 
culationscSE^ gives an exponential long-range decreasing behavior ~ G2 exp{ — |x|/T 9 }, 
where G 2 = (a s F 2 (0)) /n is the gluon condensate and T g ~ 0.2 fm Q is the gluon 
correlation length. 

In ref.,B3 the QCD two-point field strength correlator (|23|) has been related to the 
dual field propagator of the effective Abelian Higgs model describing infrared QCD. In 
this way the Gaussian dominance in the Wilson loop average is understood as following 

from the classical approximation^ in the dual theory. Moreover, it is possible to relate 
the QCD parameter T g and G2 to the dual parameters. In the London limit T g is 
identified with the dual gluon mass M, without the London limit the relation is more 
involved but still T g is expressed in terms of the dual theory parameters. 

The conclusion is that we need two parameters T g and G 2 to describe the heavy 
quark dynamics and indeed they are necessary to control the structure of the flux tube. 



Phenomenological calculation in high energy scattering indicates T g = 0.3 -7- 0.35 fm. t 
TT Surely valid in the dual description. 



Had we only one parameter, like the string tension a, we could encode the information 
of a constant energy density in the flux tube. However, the whole structure is important, 
and also the information about the width of the flux tube has to be considered. In the 
limit of very large interquark distances and in particular dynamical regimes, we can 
store the relevant information in one parameter, the string tension. 

4-3. Vacuum models for the Wilson loop 

As shown above, one obtains the 0(v 2 ) quenched quark dynamics in a given 
vacuum model, simply evaluating the Wilson loop in that model. In this way the 
phenomenological data are put in direct relation with the assumptions on the QCD 
vacuum, ilere we consider three models: Stochastic vacuum model (SVM),c3 dual QCD 
(DQCD)Ea and the flux tube model of Isgur and Paton.B^I We emphasize that these are 
models of the QCD vacuum, valid at the confinement scale^ 

• Stochastic vacuum model. The stochastic vacuum modeled is based on the idea that 
the infrared part of the QCD functional integral can be approximated by a stochastic 
process with a converging cluster expansion and a finite correlation length T g . This 
assumption is well confirmed by the lattice data. Then, the Wilson loop is given by 
Eqs.(p2|)-(p3|) with parameters T g and G2. The behavior of D and D\ and the value of 
T g are taken from the lattice measurements, G2 from the phenomenological data. The 
static potential is 

/•oo r-T 

V (r) ~ G 2 / dr d\{r - \)D(t 2 + A 2 ) (24) 
Jo Jo 

and the string tension a emerges as an integral on the D function o ~ G2 Jq 00 dr J Q °° dX 
Dir 2 + A 2 ) in the limit Tn/r — > 0. The field distribution between the quark and the 
antiquark is a fluxJsiibeo with rMs — 1-8 T 9 . Similarly the sd and vd potential are 
obtained. I refer tocj'cj for the details. Here we restrict the discussion to the long range 
limits of the potentials in the three vacuum models, (see the final subsection). 

• Dual QCD. DQCDta is a concrete realization of the Mandelstam t'Hooft dual super- 
conductor mechanism of confinement. It describes the QCD vacuum as a dual super- 
conductor on the border between type I and II. The Wilson loop approach supplies 
with a simple method to connect averaged local quantities in QCD and in dual QCD 
(=dual gluodynamics). For large loops we assumeBffil 



(W(T)) ~ (W(T)) 



Dual 
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(25) 



where C M are the dual potentials, G^ v the dual field strengths, £>j the monopole fields, 
L the effective dual Lagrangian in the presence of quarks L(G^ V ) = 2 Ti{ — \Q^ U Q^ V 
+ \{V^Bi) 2 }-U{Bi), and U{Bi) the Higgs potential. The quark sources moving around 
the Wilson loop are inserted via the Dirac string tensor G*L(x) = ge^ val 3 J ds J dty^y 13 
S(x — y). From (|25|) it follows 

gdF^izj))) = \ge, vp AG p °{z 3 ))) Dual (26) 



and it is possible to relate the averaged values of local quantities in QCD and in the 
dual theory. Then Eqs. (|15|)- dl6|) produce the 0(v 2 ) interaction in the dual formalism, 
for details cf..O We discuss only the long range limits in the final subsection. The 
nonperturbative parameters are the v.e.v. of the Higgs field, Bq, and k, the coupling 
constant of the Higgs potential (from these the penetration length and the correlation 
length can be constructed). The flux tube configurations (~ Abrikpsov-Nielsen-Olesen 
vortex) with finite tms ~ M = mass of the dual gluon, ariseEi from the numerical 
solution of the classical dual Ginzburg-Landau type of equations obtained from L. It is 
the presence of the Higgs field which confines transversally the energy distribution in a 
flux tube. Again o comes from the integration on a function exponentially decreasing 
approximately with M. 

• Isgur-Paton flux tube model. This model is extracted from the strong coupling limit 
of the QCD lattice Hamiltonian. A N-body discrete string-like model Hamiltonian 
describes the gluonic degrees of freedom. The limit N — > oo corresponds to a localized 
string with an infinite number of degrees of freedom; the radius of the flux tube is 
proportional to 1 /M. For any detail we refer to the talk of Paton in these Proceedings. 
In a recent paper,B3 the sd potential has been obtained in this model. 

4-4- Potentials and flux tube structure 

It is interesting to consider the potentials obtained using these different vacuum 
models in the limit of large interquark distance r. 
•SVM 



V (r)=ar+lcM*T g -c*T g , ± Vl {r) = -a + ^ , -fv 2 (r) = ^ , (27) 
2 ar r dr r 

V 3 and V4 fall off exponentially as exp{— r/T g } and 

. . a 1doT 2 Q 8eaT^ ta/x a doT 2 Q 2eaT„ 3 f 

Vb(r) = --r - - g - + -— ^ , K r = ~^ 9 - + -—^ , 28 

9 3 r 3 r z r 3 r z 

V d (r) = --r + -aT g -—aT g + - --JL V c (r) = --r + - -—f. 

9 4 8 3 r 9 r z b 2 r 3 r z 

• DQCD 

Vo (r) = ar-(C + C^)^- -^Vi(r) = -a + C^-~ £v 2 (r)=C~ (29) 
w v ' M dr y 1 Mr dr y ' Mr y 1 

and V 3 and V4 fall off exponentially as exp{— Mr}. 
•I-P flux tube model 

-Uf^r) = - a ^-V 2 (r) = lim ^ (30) 
dr v ' dr v ' n^oo N v ' 

and V3 and V4 are ~ 1/JV. The c, d, e, C, above are numerical constants known 
in terms of the parameters of the models. We learn that 



1) The gluon correlation length T g in SVM has the same role as the inverse of the 
dual gluon mass M in DQCD and T g /r as the inverse of the number of degrees of 
freedom N in the I-P flux tube model. All these quantities act like a correlation 
length. From now on, we use T g to refer to any of them. 



2) The MAL results are completely contained in all these mo delsft The MAL de- 

scribes the limit of large interquark distances or better the limit T g /r — > which 
is also the limit in which a potential exists, cf. Sec. 6, or the string limit. Indeed, 
all the corrections to the MAL are proportional to the gluon correlation length 
or equivalently to the width of the flux tube: the form of these corrections is 
the same in all these flux tube-like models (models that predict the flux tube 
structure) . 

3) The magnetic interaction is zero only in the limit T g /r — > 0. In the intermediate 

distance region is different from zero as one can see from the value of V2. Then a 
purely scalar effective interaction emerges only in the long range limit. 

4) T g controls the width and the shape of the flux tube as well as the validity of 

the potential description. Physically it should correspond to the size of the color 
domain as well as the fluctuation of the color fields. 

We conclude that the results for the heavy quark interaction are essentially the same 
for any flux-tube model of the QCD vacuum. 

5. Lattice calculations, flux tube models and l/Q 4 models 
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Fig. 3. Vq and V{ potentials from! 



An entire class of confining models for the heavy quark dynamics comes from the 
semirelativistic reduction of a pure convolution Bethe-Salpeter kernel (seeH and refs. 
therein) 

HQ 2 ) = iHP^UQ) + UQ) (31) 

the vd potential in DQCD see.0 
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Fig. 4. Vb and potentials from.! 
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where J v contains the one gluon exchange, P is a factor depending on the gauge and J s 
has an infrared behavior like 1/Q 4 (or 1/Q 4 , Q being the momentum transfer) chosen 
in order to reproduce the linear behavior of the static potential. To reproduce V{ then 
J s has to be scalar. These models produce: 
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VJr) 



4a a 



+ ar + C V h (r) 



3r 
V{(r) 



+ 0/— 



9r 



-(7 



y 2 '(r) 



4«s 
3r 2 



3 
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2 r T 
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(32) 



The two results in Eq. (|32] ) correspond respectively to the instantaneous form of the 
kernel 1/Q 4 or to the inclusion of retardation corrections as in£3 Notice that: 1) These 
models are constructed in order to reproduce V\ as given in Eq. (|1~8|). 2) They cannot 
reproduce the behavior of the Vy^ p (|20|) (even adding retardation corrections or allowing 
a mixture of vector and scalar confining kernelsEj). 3) The 1/Q 4 propagator does not 
avoid the spreading of the flux tube between the quarks even if the integrated energy 
increases linearly with the distance between the quarks. 

We conclude that a linear confining behavior with the correct spin-orbit is not 
sufficient to characterize the heavy quark dynamics at the order 0{v 2 ). On the other 
hand the v 2 corrections to the potential are connected to v.e.v. of field strengths inser- 
tion in the Wilson loop and therefore to the energy density and to the interquark flux 



tube structure. The 1/Q 4 models fail to take into account the dynamical nature of the 
glue. In Figs. 3-5 we present the latest lattice results for the potentials.^ The lines are 
the MAL predictions fllTf). The static potential is clearly linear for r > 0.2 fm (recent 
results indicate that the string breaking happens for r > 0.2fmcl). The nonperturbative 
behavior of V[ is consistent with the result (18) and fl32|). Unfortunately, the lattice 
data for the vd potential are not yet sufficiently accurate to discriminate between the 
flux tube model and the l/Q 4 predictions. Notice that an accurate measurement of V e 
provides a good test since it is predicted to be zero in l/Q 4 models. 

Phenomenological applications to the bottomonium and charmonium spectrum 
show that the vd corrections of the type (18) do not induce the deviations from the 
data originated by the vd corrections (32) .E3 

6. On the validity of the potential description 

In the previous sections we assumed that a potential description holds. This is 
not always the case. We know from the Lamb-shift in QED that at some order of the 
perturbative expansion the ultrasoft photons start giving a contribution to the energy 
not of a potential type. This is certainly true also in QCD but here we cannot rely 
only on the perturbative calculation. Then, the situation becomes more interesting and 
again the flux tube configuration plays a role. We can extract a potential as far as the 
adiabatic approximation holds, i.e., as far as the gluon time scale is short compared to 
the orbital period of the quarks T g «T ? . In this case we are describing a qq pair moving 
in the adiabatic potential corresponding to the ground state energy of the gluon field, 
i.e. the ground state flux tube. The next adiabatic surface describes an excited state of 
the glue: an hybrida If these adiabatic surfaces are well separated, still each of them can 
be described with a potential. The situation is controlled by the gluon correlation length 
T g which is proportional to the this difference. The quark dynamics is simply given in 
terms of the two-point correlator ( ^3|) whose nonperturbative exponential decreasing 
behavior is controlled by T g . In the case in which T g is small with respect to the other 
physical scales of the problem we recover the potential description, in the opposite limit 
we recover the sum rules description with local condensates a la Leutwyler-Voloshin.E3 
Therefore there is no discrepancy between these two descriptions. However, the cc and 
bb systems lie in the situation in which T g is small compared to the other scales. Just 
for the ground state of the bb both descriptions can be undertaken. 

Recently a new effective theory has been obtained from NRQCD integrating out 
the relative momentum scale. It is called potential Jion relativistic QCD (pNRQCD) 
and was presented at this Conference by J. Soto.EH The theory describes ultrasoft 
gluons and the matching coefficients are the potentials. pNRQCD could be a powerful 
tool to understand the relation between ultrasoft degrees of freedom and nonpotential 
nonperturbative glue. 

7. Confinement, flux tube and chiral symmetry 

As soon as the light quarks enter the game we are left without any expansion 
parameter to develop any well-founded calculation of the quark dynamics. Even in the 
simplest case, a meson formed by a light and a heavy quark, we are far away from 
the intuitive description developed for the heavy-heavy bound states. Indeed, in this 



case the heavy quark is surrounded by the fuzzy cloud of the light degrees of freedom. 
This picture is naively completely different from a flux tube picture. Moreover chiral 
symmetry breaking should appear. Here, we address the issue whether the Wilson loop 
formalism still helps to understand the quark dynamics. If this is the case we are in the 
position to study the interplay between confinement and chiral symmetry breaking. 

Let us consider the simplest system, the heavy-light bound state. Heavy-light 
systems are understood in terms of heavy quark effective theory (HQET), but HQET 
cannot predict the spectrum of the 'brown muck' since this is determined by the dynam- 
ics of strong QCD. In the limit tuq 3> A.qcd, HQET organizes a systematic expansion 
of the physical quantities in terms of Aqcd/^q and a s (rriQ). The meson mass is given 
as 

rriM = m o + A + O \ corrections. (33) 

\ m Qj 

The parameter A can be fixed on the data but its actual calculation needs a dynamical 
input. In the no- recoil limit vtlq —>■ oo, the usual choice is to take a scalar Dirac equation 
to get the A parameter as an eigenvalue. The motivation is reproducing the spin-orbit 
splitting and the fact that this kind of equation is mathematically well behaving. There 
are a number of reasons against this choice: 

• It breaks explicitly chiral symmetry. 

• We have shown that scalar confinement arises effectively for heavy quark in- 
teraction due to the collective nature of the nonperturbative gluonic degrees of 
freedom. 

In the next section we address the problem in the Wilson loop formalism. 

7.1. Gauge-invariant approach to heavy-light quark systems 

We start from the gauge-invariant quark- ant iquark Green function in the Fey- 
nman-Schwinger represent at ionE3E3: 

/•Ti m 2 _|_ ^2 

G(x,u,y,v) = -(TrP(iflp + mi ) / dT x \ Vz x e Jo 2 
4 \ y Jo Jx 

-i / dt 2 ig i dz^AJz) 

i Jo 2 e h M ; 

i rdt 1 9 -a^F^(z 1 ) i f 2 dt 2 9 -af}F^{z 2 ) . 
xe Jo 4 e Jo 4 M (_ z p v +m 2 )). 

Again the dynamics is contained in the Wilson loop, that now looks like Fig. ^. We can 
exploit the symmetry of the situation, taking the modified coordinate gauge A^(x , 0) = 
0, x^Aj(xo,yi) = (Aq(x) = Jq 1 dax h F k0 (x ,ax.), Aj(x) = Jq 1 daax k F k j(x ,ax.)) in 
which 

W(T) =TrP exp^ig dz^A^z)}. (35) 



x / dT 2 / Vz 2 e Jo 2 e Jr ~ ^ (34) 

Jo Jv 



= (-T/2,x) 
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Fig. 6. The Wilson loop in the static limit of the heavy quark and the interaction kernel K. 



At this point, at variance from the heavy quark case, we have to make a model depen- 
dent assumption: we consider still valid the dominance of the bilocal correlator with 
a finite correlation length. Under this assumption we obtain a Dirac like equation for 
the light quark in presence of the heavy quark. Approximating the Wilson loop as 



(W(T)) = ewl-££dJ*£dy»D^,y^, 
D^(x,y) = x k y l jy a a n ^ j\pp n ^\F k ,(x ,ax)FUy ,-(3y)) 



(36) 



where n(0) =0 and n(i) = 1, we obtain the Dirac-like equation for the light quark 
propagator 

S D = So + S KS D (37) 

with a kernel K(y',x') = ■y u S (y',x / )-y fl D llu (x',y'), given diagrammatically in Fig. ||. 
Notice thatEa: 1) K is not a translational invariant quantity: the coordinate gauge 
breaks explicitly this symmetry in the propagator. Physically this is due to the pres- 
ence of the heavy quark. Indeed Eq. (|37D is an integral equation for the light quark 
propagator in the field of the heavy quark. 2) The kernel depends on D^ v which in turns 
is given in terms of the two-point correlator (p3[). Then, the heavy-light dynamics is 
controlled by the same two parameters controlling the heavy-heavy dynamics, T g and 
(j2- 3) The problem has many relevant scales: the light mass m, the correlation length 
T g ~ Aqcd, the characteristic energy and momentum of the bound states. We have 
different dynamical regimes in dependences on the relative values of these scales. Let 
us study the various situations. In the following we consider only the nonperturbative 
dynamics.^ 

• Potential Case: m > 1/T g > po — m, p, p — q. We neglect the negative energy 
states and expand the kernel K in m. We obtain: 



V(r) ~ G 2 { 



dr / dX(r-X)D(r 2 +X 2 )+— — I <ir 
Jo 4m 2 r 



'2A 



r/A I — - 1J D(r 2 +X 2 )} 
(38) 



which coincides in the limit of large r with the Eichten-Feinberg potential 
with a as defined in the SVM section. We emphasize that the Lorentz structure 
which gives origin to the negative sign in front of the spin-orbit potential (hence 
to the Thomas precession term) is in our case not simply a scalar (K ~ err). 

Sum Rules case: (1/T g < po — m, 1/T g < m). We get the well-known Shifman, 
Vainshtein and Zakharov result for the heavy quark condensate 

d 4 p f d A q ro ..... 1 (aF 2 (0)) 



{QQ) = ~Ij2^I ^^WsW^oCp)} 



D s and B s case: 1/T g > m. For D s and B s one can still assume that the 
propagator inside the kernel is free and solve the equation to get the spectrum. 

D and B case: m <C l/T g . The nonlinear equation00 

S D = So + S K(S D )S D (39) 

has to be solved with Schwinger-Dyson like techniques. Notice that in the limit 
m — > the form of the kernel K is such that the interaction turns out to be 
chiral symmetric. In Ref.fcH a solution of the gap equation corresponding to Eq. 
flSUP (with a simplified form of the interaction) was found. In that case it was 
possible to disentangle the translational invariant part of the interaction (the 
self-energy) from the non- invariant part (the bound state interaction). In this 
way one obtains the light quark condensate 

(0|gg|0) ~ -T g G 2 (40) 



and the heavy-light bound state energy spectrum. The result fl40D looks appeal- 
ing: it establishes a connection between the gluon condensate and the light quark 
condensate. The connection is possible since the non-local gluon condensate has 
introduced into the game a finite correlation length T g . The same quantity that 
controls the width of the flux tube as well as the validity of the potential descrip- 
tion in the case of heavy quarks. 

In conclusion we have reached some insight in the interplay between confinement 
and chiral symmetry breaking, obtaining for the first time an unified framework for 
sum rules, potential models and chiral symmetry breaking studies. 

8. Gauge-invariant approach to light-light quark systems 



The starting point is again Eq. (pi]). The dynamics is contained in the distorted 
Wilson loop of Fig. 2. In this situation it does not exist a choice of the gauge like 
the modified coordinate gauge in the heavy-light case, i.e. it is not possible to get rid 
of both the final Schwinger strings. A result for the nonperturbative BS kernel can 
still be obtained in the form of effective diagrams^ but gauge- invariance is lost. Even 
recovering the potential from such a kernel appears to be problematic. 



9. Conclusions 



We have shown that the 0(v 2 ) heavy quark interaction is controlled only by the 
Wilson loop behavior: heavy quarks are a nice laboratory to test QCD vacuum models 
with respect to lattice and phenomenological data. On the other hand, the Wilson loop 
is a good approach to the quark dynamics: it keeps automatically the relevant degrees 
of freedom and includes naturally the flux tube configurations. The description of the 
heavy quark interaction needs two nonperturbative parameters T g and G2- these give 
origin to a only in particular dynamic situations. T g controls the width and the shape 
of the flux tube, as well as the validity of the potential description. Any model of the 
QCD vacuum reproducing the interquark flux tube structure gives the same prediction 
for the nonperturbative heavy quark interaction. On the other hand, the class of l/Q A 
models give definite different predictions. 

Moreover, the Wilson loop formalism appears to be useful in selecting the relevant 
configurations also in situations where an intuitive flux tube picture does not exist: 
e.g. the heavy-light system. In this way the heavy-heavy and heavy-light systems are 
understood in terms of the same parameters (T g , G%). These parameters are measured 
on the lattice and are also connected to the parameters of the dual superconductor 
mechanism. 

Among the open problems/work in progress we do list: extension of the analytic 
calculations to the higher order corrections in NRQCD for the heavy quark interaction; 
inclusion of the non-potential contribution (pNRQCD?); treatment in this formalism 
of systems involving more than one light quark; unquenching. 
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